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ABSTRACT 

The  plane  wave  propagation  in  a  half -space  due  to  a  uniformly  dis¬ 
tributed  step  load  of  pressure  and  shear  on  the  surface  was  first  studied 
by  Bleich  and  Nelson.  The  material  in  the  half -space  was  assumed  to  be 
elastic -ideally  plastic.  In  this  paper,  we  study  the  same  problem  for 
a  general  elastic -plastic  meterial.  The  half -space  can  be  initially 
prestressed.  The  results  can  be  extended  to  the  case  in  which  the  loads 
on  the  surface  are  not  necessarily  step  loads,  but  with  a  restricted 
relation  between  the  pressure  and  the  shear  stresses. 
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1. 


Introduction. 


Several  investigators  have  studied  the  propagation  of  elastic-plastic 
waves  of  combined  stresses  for  various  stress-strain  relations  and  for 
various  combinations  of  stress  components  [1-6].  For  general  initial 
and  boundary  value  problems,  analytical  solutions  are  difficult  to  ob¬ 
tain  except  for  a  particular  initial  and  boundary  value  problem  [2]  and/or 
a  particular  material  property  [1,6].  In  [l],  Bleich  and  Nelson  obtained 
a  closed  from  solution  for  the  case  of  plane  waves  in  an  elastic -ideally 
plastic  half-space  subject  to  a  step  load  of  pressure  and  shear  on  the 
surface.  In  [2],  Clifton  studied  the  case  of  combined  longitudinal  and 
torsional  waves  in  a  thin-walled  tube  due  to  a  step  load  of  tension  and 
torsion  at  the  end  of  the  tube.  The  materials  considered  by  Clifton  are 
general  elastic -plastic  materials  with  isotropic  work-hardening  pro¬ 

perty.  In  [6],  Ting  considered  the  case  of  two  shear  waves  in  an 
elastic,  linearly  work-hardening  half -space  subject  to  a  series  of  step 
loadings  and  unloadings • 

In  the  analyses  of  wave  propagation  in  a  thin-walled  tube,  the 
lateral  inertia  was  ignored.  While  this  is  a  good  approximation,  care 
must  be  exercised  in  comparing  the  theoretical  analyses  with  an  experi¬ 
mental  result.  For  wave  propogation  in  a  half -space,  the  problem  of 
lateral  inertia  does  not  arise.  Thus  the  results  of  pressure -shear 
waves  studied  by  Bleich  and  Nelson  can  be  used  without  reservations  for 
experimental  verifications.  Since  an  elastic -ideally  plastic  material 
is  an  idealization  of  real  materials,  the  study  of  pressure-shears  waves 
in  a  half -space  of  general  elastic -plastic  materials  is  desirable. 

This  is  presented  in  this  paper. 
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The  governing  differential  equations  of  the  problem,  the  character¬ 
istics  and  the  eigenvectors  are  presented  in  section  2-  In  section  3> 
we  obtain  the  stress  paths  in  the  stress  space  which  are  used  to  obtain 
simple  wave  solutions.  It  is  shown  by  an  example  that  a  simple  wave 
solution  can  be  a  solution  for  the  case  in  which  the  pressure  and  the 
shear  applied  at  the  boundaiy  are  arbitrary  with  a  restricted  relation 
between  the  pressure  and  the  shear.  The  particular  case  in  which  the 
material  is  elastic -ideally  plastic  as  considered  in  [1]  is  reduced  in 
section  h,  but  presented  in  a  form  which  can  be  used  also  when  the  half- 
space  is  initially  prestressed. 


?.  The  Basic  Equations. 

Let  the  half-space  be  bounded  by  the  horizontal  plane  x  *  0  of 
tm*  cartesian  coordinates  (x,y, z)  and  extended  to  infinity  on  the  side 
for  which  x  >  0.  Let  u(x,t)  and  v(x,t)  be  the  x  and  y  components 
of  the  velocity  of  any  particles  which  depend  only  on  x  and  the  time 
t.  The  z  component  of  the  velocity  is  assumed  to  be  zero.  The  equa¬ 
tions  of  motion  for  this  plane  motion  is 


(1) 


dv 

cfc  =  p  dt 


(2) 


where  p  is  the  mass  density  of  the  half-space,  and  a,  -  o  ,  i  =  t 

j.  xx  xy 

for  simplicity.  The  stress-strain  relation  for  an  isotropic  work-harden¬ 
ing  material  is  (see  (73) 


where  E  is  Young's  modulus*  v  is  Poisson's  ratio  and  \  is 
a  parameter  which  will  be  defined  shortly,  f  is  the  yield  condition 
which  can  be  written  as 


f  =  ( 


Vff2-2 


)  +  T2  =  k2 


(4) 


where  =  a  =  a  *  k  is  the  yield  stress*  and  9  is  a  constant 
2  yy  zz 

which  assumes  the  value  -J5  for  the  von  Mises  yield  condition  and  the 
value  2  for  the  Tresca  yield  condition.  Applying  Eqs.  (3)  and  (4)  to 
the  problem  under  consideration*  and  noticing  that  =  0  and 


-v  cie 
ou  xx 


_  1  dv  d£Xy 

sr~5T'  2  Sc  =  5t 


(5) 


by  the  continuity  requirement*  we  obtain 

du  1  S(J1  2v  Sa2  b\ 
ST  I3t - e3T  sSt 


(6) 


0  = 


2v  S(I1  ,  2(l-v)  ^a2 

1F+  e  S'  “  s  at 


(7) 


Sv 

Sc 


1  dx 
p  St 


-I 


where  p  is  the  shear  modulus  and 


s 


02  (ar02}. 


(8) 


(9) 


X  in  Eq.(3)  can  be  expressed  in  terms  of  f  as  (see  [7])* 


(10) 


8x 

St  = 


e2o(k)  5f  s“it i 


where  a(k)  characterizes  the  work-hardening  property  which  can  be 

determined  from  the  stress-strain  curve  for  a  simple  tension  test.  If 

E  (cr)  is  the  slope  of  the  stress-strain  curve  for  a  simple  tension  test 
P 

expressed  in  terms  of  the  tensile  stress  a ,  we  have  (see  [2]), 


Ep(0k)  "  1  ^i:L^ 

In  the  elastic  region,  E  =  E  and  a  =  0,  while  in  the  ideally  plastic 

P 

region,  E^  =  0  and  a  =  «.  We  shall  assume  that  E^  is  a  monotonically 
decreasing  function  of  k.  Then,  by  Eq.(ll),  a(k)  is  an  increasing 
function  of  k. 

Equations  (l),  (2),  (6),  (7),  (8),  and  (10)  can  be  written  in  a 
matrix  form: 


A  w,  +  B  w  =0 
~  ~t  ~  ~x 


(12) 


where  the  subscripts  x  and  t  denote  partial  differentiation  with 
respect  to  these  variables  and 


P  0 
0  p 
0  0 
0  0 
0  0 
0  0 


0  0 

0  0 

1  _  2v 

E  E 

2v  2(l-v) 
E  E 

0  0 

s  -s 


0 

0 

0 

0 

1 

n 

2t 


0 

0 

s 

-s 

2t 

4k2E 

e2a 
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0  0  -1  0  0  0 


0  0  0  0  -1  0 


-1  0  0  0  0  0 


0  0  0  0  0  0 


0  -1  0  0  0  0 


0  0  0  0  0  0 


Notice  that  both  matrices  A  and  B  are  symmetric .  By  letting  a  -» «, 
we  obtain  alternate  squat ions  for  elastic -ideally  plastic  materials 
considered  in  [ l] .  The  characteristics  c  of  Eq.(l2)  are  the  roots  of 
the  determinant  |cA-b|  =  0 ,  (see  [8]).  After  expanding  the  determinant. 


we  obtain: 


c  D(c)  =  0 


where 


d(c)  =  ~ 

2  p  p  92a(p+;_) 


(14a) 


-(f-)2((l4)s2  +  ^  t2+%^J+(s2+%^)} 
2  P  0  2Q(p+l)  02Ct(  p+l) 


2 

Cg  =  |j/p  is  the  shear  wave  speed  and  P  is  the  parameter  introduced  by 
Bleich  and  Nelson  which  is  related  to  Poisson's  ratio  v  by  the 


equation 


P=2ilivi 

p  1-2V  k  ** 

As  v  assumes  the  range  (0,  p  has  the  range  (2,oo).  From  Eq. 

(15),  we  have  either  c  =  0  or  c  obtained  by  the  roots  of  D(c)  =  0. 
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c„  and  c  correspond,  respectively,  to  the  fast  wave  speed  and  the 
I  s 

slow  wave  speed. 

Two  ex+reme  cases  can  be  reduced  from  Eq.(l4b).  In  the  elastic 
region,  a  =  0  and  D(c)  =  0  gives  c  =  c^  and  c  =  c-^.  In  the 
ideally  plastic  region,  O.  =  »  and  the  roots  of  Eq.(l4b)  reduce  to  the 
ones  obtained  in  [ l] . 

Other  extreme  cases  can  be  reduced  from  Eq.(l4a).  We  simply 
list  the  results  in  the  following: 

(i)  When  s  =  0, 

°f  *  °1  '  0  ?.  %  -  °2  (l8> 
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(ii)  When  t  =  0,  and  0  >  3, 


then  Eqs.  (4)  and  (9)  give 

p 

s  =  q  k  cos  <p  .  (22b) 

With  this  change  of  variables,  cf  and  cg  become  functions  of  k 

and  <p.  In  the  s  ~  'r  plane,  a  constant  k  gives  a  yield  surface  while 

a  constant  <p  gives  a  straignt  line  through  the  origin.  Regarding  cf 

and  c  as  functions  of  k  and  <p,  it  can  be  shown  that,  for  0  <  q>  < 
s 


8 


ro|a 


8c 

8c 

8c 

8c 

w>0 ' 

W<0  > 

5F<o- 

sf<° 

(23) 


From  Eq.(23),>  and.  the  information  presented  in  Eqs.(l8)  to  (20), 

and  c  have  the  ranges: 
s 


VI 

o 

°s  S  C2  ’ 

I  °2  —  Cf  —  C1 

for 

P  >  3 

VI 

o 

%  -  C2  ' 

C2-Cf-  C1 

for 

P  <  3 

c 


f 


(24a) 

(24b) 


The  characteristic  condition  along  a  characteristics  c  is  obtained 


by  [8], 


„T. 


dw 


i!  A  -nr  =  0 
~  ~  dt 


(25) 


r*J 

where  —  =  is  the  total  derivative  along  the  characteristics 

m 

c  and  i-  is  the  transpose  of  the  left  eigenvector  £  which  is  obtained 
by  the  equation 


(cA-B)  =  0  . 

'  rst* 


(26a) 


In  particular,  for  the  characteristics  c  =  0,  Eq.(25)  reduces  to  Eqs. 
(7)  and  (10). 

When  A  and  B  are  symmetric  as  in  the  present  case,  Eq.(26a)  can 
be  written  as 


(cA-B)£  =  0.  (26b) 

Thus  the  left  eigenvector  and  the  right  eigenvector  are  identical. 

For  c  =  +  c„  or  +  c  ,  £  of  Eq.(26b)  has  the  solution 
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r 


where 


! 

i 

&  =  I  -PC$ 
<5 

-pc 

P> 


*  _  l!l  .  sP  n  .6*4  4k2..w-  _% 

kT  ^  p+i  a2  2 

K  0  s  0!  C 


s  C2 
c2 


* =  lfc§l  £si  +  <*§§.  (  +  ti.  JaL)(1.! I) 

3  s  4t  U  p+le2  2aAA  c2' 


(27) 


(28a) 

(28b) 


°2 


$  and  4>  can  be  writeen  in  simpler  forms  by  eliminating  a  from  the 
equation  D(c)  =  0  where  D(c)  is  expressed  in  Eq.(l4b): 


2  2 
s  c  "c2 

^  =  t  ~2  2 
c  -cx 


<J>  = 


Ps(c2-c2)(c2-Pc2) 

2tc(c2-c2) 


(29a) 


(29b) 


3*  Generalized  Simple  Wave  Solutions. 

Simple  wave  solutions  are  particular  solutions  of  Eq.(12)  in  which 
w  is  a  constant  vector  along  a  characteristics  [9] •  Thus,  if  w  is 
a  constant  vector  along  the  c^  we  have  a  fast  simple  wave  solution. 
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If  w  is  a  constant  along  the  c  we  have  a  slow  simple  wave  solution 

s 

[2] .  According  to  the  theory  of  generalized  simple  waves  developed  in 
[9];  in  the  region  where  the  simple  wave  solution  is  valid,  dw  is 
proportional  to  the  eigenvector  £.  Therefore,  by  Eq.(27), 


du 

$ 


dv 

1 


d°l 

-pci 


dar 

C 

0 


dr 

-pc 


dk 

© 


(30) 


Equation  (30)  is  equivalent  to  five  differential  equations  of  the  first 
order.  Since  J,  $  and  ©  are  functions  of  the  stresses  only,  these 
five  equations  are  not  all  coupled.  In  particular,  since 


da  i“d<Tg 


dr 


-pc$-<t  -pc  ’ 


we  have,  by  Eqs.(9)  and  (28), 


12pk 


dS 

dT  e2s  e4TS(p+i)a 


(1-4)  , 


(31a) 


or,  by  Eqs.(9)  and  (29), 


as  s(c2-=|)(3=a-Pc|) 

dT  fl2  2/2  2n 
9  TC  (c  -c^) 


(31b) 


The  right  hand  side  of  Eq.(31a)  or  Eq.(31b)  is  a  function  of  s  and 
t  only.  Hence  Eq.(31a),  or  Eq.(31b),  is  a  differential  equation  by 
itself.  To  obtain  the  stress  paths  in  the  (a^, ct^,t)  space  for  simple 
wave  solutloi. „,  we  need,  in  addition  to  Eq.(51b),  another  equation  which 
can  be  written,  by  Eqs.(30)  and  (29a), 


da.l 

dr 


2  2 

s  C  ~C2 

T  2  2  ' 

c  -C]_ 
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(32) 


Equations  (31b)  and  (32),  when  integrated,  given  two -parameter  family 
of  space  curves (for  each  o,  and  c  )  in  the  three-dimensional  space 

A  S 

(ai,Og,T).  These  space  curves  can  be  constructed  if  we  have  their 
projections  on  any  two  planes  which  are  not  parallel  such  as  the 
and  the  t  planes.  In  view  of  Eqs. (31b)  and  (32),  we  will  use  the 
s~t  plane  and  the  a^~  t  plane  instead.  The  reason  is  obvious. 

Equation  (31b),  when  integrated,  gives  only  one-parameter  family  of 
curves  in  the  s~t  plane.  As  to  the  projected  curves  on  the  a. '-'T 
plane,  we  also  need  only  one -parameter  family  of  curves.  The  second 
parameter  family  of  curves  can  be  obtained  by  translating  the  first 
parameter  family  of  curves  in  the  (indirection.  This  is  clear  from 
Eqs •  (31b)  and  (32). 

Before  we  illustrate  how  to  construct  these  curves  (or  stress  paths), 
for  simple  wave  solutions,  we  will  present  some  properties  of  these 
curves  in  the  s~r  plane  and  the  o1~T  Plane-  First,  consider  these 
curves  in  the  a^-^  plane.  By  Eqs.  (32)  and  (17),  we  have 


da 

(~)  <  0 
dr  /c=cf 


da. 

<dTU  >0 

s 


(33) 


Using  Eqs.  (32)  and  (l4b),  it  can  be  shown  that 

do.  da.. 

( — i)  .  ( — i)  = 

'dT  'c-C,  '/,'r  'r'~r 


'dr  /c=c 


m 


f  s 

In  other  words,  the  stress  paths  for  the  fast  simple  waves  and  the  slow 
simple  waves  are  orthogonal  to  each  other  in  the  o^  ~  t  plane.  It 
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should  be  noticed  that  they  are  orthogonal  in  the  ~  t  plane  provided 
they  intersect  each  other  in  the  (a^,  ct^t)  space.  If  a  stress  path 
for  the  fast  simple  wave  and  a  stress  path  for  the  slow  simple  wave  do 
not  intersect  in  the  space,  their  projections  in  the  ~  r 

pxane  are  not  necessarily  orthogonal  to  each  other. 

Next,  consider  the  projections  of  the  stress  paths  in  the  s~t 
plane.  From  Eqs.  (24)  and  (31a),  we  have 


(^)c=c  ~T  ’  for  a11  p  • 

a  c  cf  @  S 

2/  2 


(35) 


2  2 

For  slow  waves,  ds/dr  depends  on  whether  CJC2  ^ar6er  or  smaller 

2  2  / 
than  p/3  •  (if  P  >  3,  c/c  is  automatically  smaller  than  P/3  • )  By 

S  ed 

Eqs-  (24),  (31a),  and  (31b),  we  have 


4t  ./ds\  ^  _  _  Cs  ^  3 

fl2  -  ^dx^c=c  -  0  '  for  2  <  3  * 

os  s  c2 


(36a) 


s  8 

If  -p  >  £,  (this  can  happen  only  when  P  <  3)  it.  is  seen  by  Eqs. 
c  ^ 


(24b)  and  (31b)  that 


0  <  (~)  <  -  ,  for  —  >  | 

-  dT'c=C  -  T  ’  2-3 

s  c„ 


(36b) 


In  Figs-  1  and  2,  we  give  two  examples  of  stress  paths  projections 
in  the  ~  t  plane  and  the  s  ~  t  plane  by  integrating  Eqs.  (3lb) 
and  (32)  numerically.  s  in  Figs.  1(b)  and  2(b)  has  been  changed  to 
(ci~a2)  according  to  Eq.  (9)«  kQ  is  the  initial  yield  stress.  In 
Fig.  1,  p  =  5,  which  shows  a  typical  example  of  p  >  3  while  in  Fig.  2, 
the  value  p  =  2.25  gives  an  example  of  P  <  3*  In  both  cases,  the 
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2 

von  Mises  criterion  is  used  and  hence  9  =  3*  The  work-hardening  pro¬ 

perty  a(k)  are  chosen  differently  in  Figs.  1  and  2.  They  are  chosen 
for  illustrative  purpose  and  do  not  necessarily  represent  a  real  material. 
In  both  figures,  the  solid  lines  marked  by  s^, s^, . ..,  are  the  stress 
paths  for  the  slow  simple  waves  while  the  dashed  lines  marked  by 
f 1* fg, •  •  * ,  are  the  stress  paths  for  the  fast  simple  waves.  The  sub¬ 
scripts  1,2,...,  have  no  particular  meaning  except  to  identify  the 
stress  paths  in  the  (o^, o^,?)  space.  For  Instance,  the  curve  s^  in 
Fig.  1(a)  and  the  curve  s^  in  Fig.  1(b)  represent  the  same  space 
curve  in  the  t)  space.  In  Fig.  2(b)  the  point  d  on  the  s-axis 

corresponds  to  the  value  s*  defined  in  Eq.  (2l).  In  fact,  s  =  s*, 
t  =  0  is  a  singularity  of  E'qs.  (31b)  and  (32).  The  dotted  line  in 
Fig.  2(b)  is  the  locus  of  c2/c2  =  p/3«  It  is  seen  that  Eqs.  (33)~(36) 

S  cL 

are  satisfied  by  the  results  obtained  in  Figs.  1  and  2. 

The  arrow  heads  in  the  figures  show  the  directions  along  which  the 

stress  state  should  be  changed  to  insure  a  continuous  plastic  flow. 

It  is  apparent,  by  Eq.  (23)  and  Eqs.  (33) -(36),  that  along  these 

directions  shown  by  the  arrow  heads,  the  wave  speeds  c„  and  c 

±  s 

decrease  on  the  strese  paths  for  the  fast  simple  waves  and  the  slow 
simple  waves  respectively.  This  is  an  important  requirement  for  con¬ 
structing  a  simple  wave  solution. 

In  integrating  the  differential  equations,  Eqs.  (31b)  and  (32) 
initial  conditions  are  required.  The  initial  conditions  for  the  results 
obtained  in  Figs.  1  and  2  are  chosen  in  such  a  way  that  Figs.  1  and  2 
can  be  used  to  construct  simple  wave  solutions  for  the  case  in  which  the 
half-space  is  initially  stress  free.  No  new  calculations  are  required 
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if  the  half-space  is  initially  pre-stressed.  By  rearranging  the  curves 
in  Figs.  1(b)  and  2(b),  Figs.  1  and  2  can  be  used  for  the  case  in  which 
the  half-space  is  initially  pre-stressea-  We  will  illustrate  this  later 
in  this  section. 

When  the  half-space  is  initially  stress  free,  the  disturbance 
caused  by  a  load  on  the  surface  of  the  half -space  will  propagate  at  the 
elastic  wave  speed  c^  or  c^*  The  response  at  the  wave  front,  and 
possibly  a  finite  region  behind  the  wave  front,  is  elastic  and  Eqs. 

(1),  (2),  (6),  (7),  and  (8)  apply  if  we  let  dx/dt  =  0.  In  particular, 
Eq.  (7)  yields 


v 

1-v 


or,  since  cr,  and  cr_  are  initially  zero,  we  have,  making  use  of 
1  d 

Eq.  (15), 


_  P-2 

°2  p+4  al  * 


(37) 


Therefore,  the  initial  yield  limit  kQ  is  reached  when,  by  Eq.  (4), 


<^)2  +  T*  =  £  ■ 


(38) 


Equation  (58)  is  the  initial  yield  surface  shown  in  Figs.  1(b)  and  2(b) 
while  elimination  of  between  Eqs.  (57)  and  (58)  gives  the  initial 
surface  shown  in  Figs.  i(a)  and  2(a).  In  other  words,  Eqs.  (57)  and 
(58)  represent  a  curve  in  the  space  and  the  projections  of 

this  curve  in  the  s  ~  t  plane  and  the  a1  ~  t  plane  are  denoted  by 
initial  yield  surface  in  Figs.  1  and  2.  The  initial  conditions,  or  the 
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"starting  points",  for  the  stress  paths  for  the  fast  simple  waves  shown 
in  Figs.  1  and  2  are  taken  on  the  T-axis  and  on  the  curve  represented 
by  Eqs.  (37)  and  (38).  The  starting  points  of  the  stress  paths  for  the 
slow  simple  waves  are  taken  on  the  curve  represented  by  Eqs.  (37)  and 
(38),  and  on  a  curve  on  the  a ^  ~  a g  plane  in  the  (a^a^T)  space 
which  we  will  derive  next.  This  curve  is  in  fact  a  stress  path  for 
the  fast  simple  waves  which  lies  on  the  ~  a ^  plane. 

From  Eqs.  (31b)  and  (32)  we  have,  for  c  =  c„, 


ds  _  3  P  °f 
dal  6 2  e2  c2 


(39) 


0  2  2 

When  t  =  0,  —  s  =  k  and,  after  solving  cp/c  from  Eq.  (13)  and 

2  x  c. 

substituting  the  result  into  Eq.  (39);  we  obtain 


0  ds  _ _ 6 _ _  _  /i, 

2  da1  =  (P+l)a(0s/2)+(P+4)  for  T  =  0  ' 

This  is  the  differential  equation  for  the  stress  paths  for  the  fast 
simple  waves  when  t  =  0.  Since  a(k)  is  a  given  function  of  k, 
a(e  s/2)  is  a  known  function  of  s  and  Eq.  (40a)  can  be  integrated. 
The  initial  condition  is  taken  from  the  intersection  of  the  curve 
represented  by  Eqs.  (37)  and  (38)  with  the  plane  t  =  0,  which  gives: 


s  =  |  ;  &i  =  |  (P+4)kQ  .  (40b) 

Equations  (40a)  and  (40b)  represent  a  curve  on  the  plane  i  =  0.  The 
starting  points  of  the  stress  paths  for  the  slow  simple  waves  shown  in 
Figs.  1  and  2  are  taken  on  the  curve  represented  by  Eqs.  (57)  and  (38) 
and  the  curve  represented  by  Eqs.  (40a)  and  (4ob). 
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Now,  we  are  in  the  position  to  construct  simple  wave  solutions  by- 
using  the  stress  paths  obtained  in  Figs.  1  and  2  for  the  case  when  the 
half -space  is  initially  stress  free.  The  stress  paths  shown  in  Figs.  1 
and  2  are  for  a  plastic  region.  In  an  elastic  region,  the  stress  paths 
corresponding  to  the  fast  simple  waves  c  =  c1  are  horizontal  lines 
parallel  to  tne  and  &v&xes  while  the  stress  paths  corresponding 
to  the  slow  simple  waves  c  =  are  vertical  lines  parallel  to  the 
T'-axis.  The  stress  paths  for  an  elastic  region  are  not  oriented,  i.e., 
the  stress  state  can  be  changed  in  either  direction.  We  define  an 
"admissible  stress  path"  by  the  path  which  consists  of  one  or  more 
stress  paths  in  the  elastic  and/or  plastic  regions  in  such  a  way  that 
when  one  moves  along  the  path,  the  wave  speeds  are  non-. ncr easing. 

Thus,  for  instance,  the  paths  oaij,  obg  and  ode  in  E.'g.  2  are 
admissible  stress  paths.  For  each  admissible  stress  path,  a  simple 
wave  solution  can  be  constructed.  For  instance,  a  simple  wave  solution 
for  the  admissible  stress  path  obg  is  shown  in  Fig.  3(a)  where  the 
solution  corresponds  to  that  of  a  step  load  of  cr^  =  o®,  t  =  on  the 
surface  of  the  half-space.  A  superscript  g  denote  the  value  at  the 
point  g  in  Fig.  2.  Since  the  characteristics  c  are  functions  of 
the  stresses  cr^,  a^,  t,  and  since  the  stresses  are  constant  along 
the  characteristics  for  a  simple  wave  solution,  the  characteristics  are 
straight  lines  as  shown  in  Fig.  3(a).  Each  point  on  the  stress  path 
obg  corresponds  to  a  characteristic  line  in  the  x  ~  t  plane.  The 
position  of  the  point  on  obg  determines  the  stresses  along  the 
characteristics  and  also  the  slope  of  the  characteristics.  Thus 
and  t  are  determined  from  Fig.  2(a)  and  is  determined  from  Fig. 
2(b)  accordingly. 
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Fig.  3(b)  shows  another  simple  wave  solution  corresponding  to  the 
same  stress  path  obg.  The  loads  on  the  surface,  a^(0,t)  and  T(0,t) 
are  not  step  loads.  cr^(0,t)  and  T(0,t)  however,  should  be  prescribed 
in  such  a  way  that,  as  t  increases,  cr^(0,t)  and  f(o,t)  follow  the 
stress  patn  obg.  Thus,  one  of  them,  say  a^(0,t),  can  be  prescribed 
almost  arbitrarily,  and  the  other  one  T(0,t)  is  obtained  by  cr^(0,t) 
and  the  stress  path  obg.  It  should  be  noted  that  both  a^(0,t)  and 
T(o,t)  do  not  have  to  be  continuous  functions  of  t  as  shown  in  Fig. 
3(b). 

When  the  half-space  is  prestressed,  the  stress  paths  in  the  s  ~  t 
plane  remain  unchanged  while  the  stress  paths  in  the  cr  ~  t  plane  need 
some  modifications.  To  illustrate  this,  let  us  consider  the  case  in 
which  the  half -space  is  pre-sheared.  In  Fig.  1,  suppose  that  the  pre¬ 
sheared  stress  is  at  the  point  a.  An  admissible  stress  path  in  Fig. 
l(b)  is  abdegh  which  consists  of  the  stress  paths  f0  and  s^.  The 
corresponding  stress  path  in  Fig.  l(a)  however,  is  not  ab'd' e' g'h' 
since  g  in  Fig.  l(b)  and  g'  in  Fig.  l(a)  do  not  represent  the  same 
point  in  the  (a^a^r)  space.  This  can  be  verified  easily  as  the 
ordinates  of  g'  and  g  are  not  the  same.  Similarly,  b' ,  d' ,  e' 
in  Fig.  l(a)  and  b,  d,  e  in  Fig.  1(b)  do  not  represent  the  same 
points  in  the  (a^,Og,T)  space.  As  we  explained  earlier,  the  stress 
paths  obtained  in  the  <7j~t  plane  can  be  used  for  other  initial  condi¬ 
tions  if  we  translate  the  curves  in  the  horizontal  direction.  Thus  we 
translate  the  curves  s^  Sg,  s^,  and  s ^  in  Fig.  1(a)  until  their 
intersections  with  fg  give  the  same  ordinates  as  b,  d,  e  and  g 
respectively  of  Fig.  l(b).  The  result  of  this  translation  is  shown  in 
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Fig.  4.  Now,  abdegh  in  Fig.  l(b)  and  Fig.  4  represent  a  continuous 
curve  in  the  (o^o^x)  space,  so  do  abdi  and  abdegj .  They  are  all 
admissible  stress  paths.  For  the  stress  path  abdegj,  the  last  segment 
gj  corresponds  to  an  elastic  unloading.  Once  we  obtain  an  admissible 
stress  path,  a  simple  wave  solution  can  be  constructed  as  illustrated 
in  the  previous  example,  and  Figs .  3  * 

The  case  when  the  half-space  is  pre-compressed  can  be  analysed  in 
a  similar  manner.  If  the  half -space  is  pre-sheared  and  pre-compressed, 
then,  in  addition  to  the  initial  values  of  and  t,  the  initial 
value  of  a2  has  to  be  specified.  This  is  so  because  depends  on 
how  the  half-space  is  pre-sheared  and  pre-compressed.  Again,  Figs.  1 
and  2  can  be  used  for  this  case  with  minor  modifications. 

4.  Ideal  Elastic -Plastic  Materials. 

When  the  material  is  elastic -ideally  plastic,  i-e.,  a  =  “,  the 
analysis  is  greatly  simplified.  Instead  of  projecting  stress  paths  on 
the  s  ~  t  plane  and  the  0.^  ~  t  plane,  only  the  latter  is  required. 
Moreover,  only  one  stress  path  for  the  fast  simple  waves  and  the  slow 
simple  waves  need  to  be  calculated;  the  rest  are  obtained  by  a  transla¬ 
tion  of  the  curve  calculated.  No  modifications  are  necessary  for  the 
case  when  the  half-space  is  pre-streased. 

For  an  ideally  elastic-plastic  material,  the  yield  condition  Eq.  (4) 
is  replaced  by  Eq.  (38)  where  kQ  is  a  fixed  constant.  Thus  0g  is 
no  longer  an  independent  unknown  but  can  be  expressed  in  terms  of  0^ 
and  t.  The  characteristic  equation  D(c)  =  0  now  becomes,  by  Eqs. 

(l4b)  and  (38), 
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,.2  2 w  2  2U,  2  a  2,  .  .2  2  2,  2  2,  . 

(kQ-T  )(c  -C2)(3c  -PCg)  +  6  T  C  (C  "C^)  =  0 


(4l) 


Therefore,  c  does  not  depend  on  The  values  of  c  when  t  =  0 

and  t  =  kQ  can  be  obtained  easily  from  Eq.  (4i).  The  result  is 


Cf  =  JJ  c2  *  cs  =  c2  >  when  T  "  0  and  P  >  3 


(42a) 


Cf  =  0o  ,  Cg  =  I-  cg  ,  when  t  =  0  and  p  <  3 


(42b) 


c.  =  c..  ,  c  =  0  ,  when  t  =  k  . 
X  X  s  o 


(42c) 


It  can  be  shown  that 


Sc  Sc 

%r<0  >  ^r>  o 


(43) 


Hence,  by  Eqs.  (42)  and  (43),  c„  and  c  nave  the  ranges 


0  <  Cs  -  C2  *  J3  °2  -  °f  -  C1  for  P  -  3 


(44a) 


0  <  c 
—  s 


<  Cg  >  c2  S  cf  —  ci  f°r  P  <  3 


(44b) 


This  agrees  with  the  result  obtained  in  [l]. 

Using  Eqs.  (38)  and  (4l),  the  right  hand  side  of  Eq.  (32)  can  be 
written  in  terms  of  a,  and  t  alone.  The  result  is, 

i  * 


da1  Q  +  iQ2+l6©2r2(k^-T2) 


dT 


40  t  /k2--^ 
1  o 


where 
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Q  =  0~3)k2  +  02  -(p-3)) t2 


The  +  sign  in  Eq.  (45)  is  for  the  fast  simple  waves  while  the  -  sign 
is  for  the  slow  simple  waves.  Equation  (45)  can  be  integrated  in  a 
closed  form  in  terms  of  elliptic  integrals-  The  particular  case  in 
which  0  =3  was  obtained  in  [l].  In  Fig*  5>  we  show  the  results  of 
an  integration  of  Eq.  (45)  for  three  cases:  P  >  3>  P  =  3>  and  P  <  3> 

In  each  case  only  one  curve  for  the  fast  simple  waves  and  the  slow 
simple  waves  needs  to  be  calculated.  The  rest  of  the  curves  are  obtained 
by  translation  in  the  (^direction  of  the  curve  calculated.  Again, 
the  solid  lines  with  arrow  heads  are  the  stress  paths  for  the  slow 
simple  waves  while  the  dashed  lines  with  arrow  heads  are  for  the  fast 
simple  waves.  In  addition,  the  o^axis  is  also  an  admissible  stress 
path  which  corresponds  to  a  constant  wave  speed  of  J p/3  c^.  Notice 
that  Eq.  (34)  still  applies  and  therefore  the  stress  paths  for  the 
fast  simple  waves  and  the  slow  simple  waves  are  orthogonal  to  each 
other  everywhere*  Notice  also  the  angle  at  which  the  stress  paths 
intersect  the  a^axis.  The  stress  paths  for  the  slow  simple  waves 
shown  by  solid  lines  intersect  the  a^~axis  at  90°  for  P  >  5,  at  45° 
for  P  =  3,  but  never  intersect  the  a^~axis  for  p  <  3-  On  the  other 
hand,  the  stress  paths  for  the  fast  simple  waves  shown  by  dashed  lines 
never  intersect  the  o^~axis  for  P  >  5  but  intersect  the  a^~axis  at 
45°  for  P  =  3  and  at  90°  for  P  <  3*  The  stress  paths  for  fast  and 
slow  simple  waves  do  inteiseet  the  horizontal  line  t  -  kQ>  one  inter¬ 
sects  at  90°  and  the  other  is  tangent  to  the  line  t  =  k  . 
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Fig-  5  can  be  used  for  obtaining  simple  wave  solutions  regardless 
of  whether  the  half -space  is  initially  stress  free  or  not-  Any  point  in 
the  ~  t  plane  can  be  taken  as  the  initially  prestressed  state  aid 
an  admissible  stress  path  is  then  determined  accordingly.  For  instance, 
if  the  point  a  is  the  initially  prestressed  state,  the  path  abde  is 
an  admissible  stress  path.  So  is  the  path  agh  for  (3  >  3  and  (3  -  3 
but  there  is  no  corresponding  path  for  P  <  3  since  the  solid  lines 
never  intersect  the  cr^axis.  If  point  d  in  Fig.  5(a)  is  the  initially 
prestressed  state,  the  stress  paths  dih  and  djm  are  both  admissible. 

As  a  last  example,  suppose  that  the  half -space  is  initially  stressed 
beyond  the  initial  yield  surface  as  shown  by  the  point  i  in  Fig.  5(c), 
and  at  t  =  0  a  step  load  of  is  added  on  the  surface  of  the  half- 

space  such  that  the  stress  state  changes  to  the  point  j  in  Fig.  5(c). 

The  admissible  stress  path  for  this  case  is  igj  which  shows  that  the 
shear  stress  will  decrease  first  and  increase  later  even  though  the 
shear  stress  on  the  surface  of  the  half-space  is  kept  constant. 

Once  we  have  an  admissible  stress  path,  a  simple  wave  solution  can 
be  constructed  following  the  examples  shown  in  Figs.  3- 

5 •  Concluding  Remarks. 

The  problem  of  plane  waves  of  combined  pressure  and  shear  stresses 
in  an  ideally  elastic -plastic  half -space  originally  studied  by  Bleich 
and  Nelson  is  extended  to  a  general  elastic -plastic  material.  Simple 
wave  solutions  are  obtained  for  the  half-space  which  can  be  initially 
pie-stressed  and  the  loads  on  the  boundary  need  not  be  a  step  load. 

The  simple  wave  solutions  presented  here  of  course,  do  not  apply  to 


22 


general  boundary  value  problems.  For  general  boundary  value  problems, 
a  numerical  scheme  such  as  the  method  of  characteristics  using  Eq.  (25) 
must  be  used.  The  simple  wave  solutions  obtained  here,  such  as  the  one 
shown  in  Fig.  3(b),  can  be  used  as  a  test  of  the  accuracy  of  the  numeii- 
cal  scheme. 

The  existence  of  the  fast  waves  and  the  slow  waves  has  been  verified 

experimentally  by  Lipkin  and  Clifton  [10].  An  example  in  which  all  the 

four  wave  speeds  cn,  c_,  c„  and  c  are  generated  dn  one  test  is 

l  d  f  s 

given  in  Fig.  6.  For  illustrative  purpose,  we  assume  the  half -space  is 

elastic -ideally  plastic.  The  half -space  is  initially  pre-stressed  to 

the  stress  state  indicated  by  point  a  in  Fig.  6(a),  and  at  time  t  =  0 

the  stress  stat;  at  the  boundary  is  changed  to  the  stress  sate  indicated 

by  point  g  in  Fig.  6(a)  and  maintains  at  this  constant  state  thereafter. 

For  this  case,  the  admissible  stress  path  is  abdeg.  The  portion  ab 

is  elastic,  with  wave  speed  c^j  the  portion  bd  is  plastic,  with 

variable  wave  speeds  c^;  the  portion  de  is  again  elastic,  with  wave 

speed  Cgj  and  the  last  portion  eg  is  plastic,  with  variable  wave 

speeds  c  .  The  corresponding  simple  wave  solution  is  shown  in  Fig. 
s 

6(b). 

In  the  wave  propagation  of  combined  stresses,  an  unexpected  un¬ 
loading  may  occur  near  the  boundary  when  the  stress  state  at  the  boundary 
suddenly  changes  from  a  lower  yield  surface  to  a  higher  yield  surface, 

[2] .  On  the  other  hand,  a  plastic  loading  may  occur  when  the  stress 
state  at  the  boundary  suddenly  changes  from  a  higher  yield  stress  to  a 
lower  yield  stress,  such  as  the  path  djm  in  Fig.  5(a).  This 
phenomenon  is  not  dictated  by  the  fact  that  the  stress  state  at  the 
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boundary  changes  discontinuously .  The  same  phenomenon  exists  if  the 
stress  state  at  the  boundary  changes  continuously.  In  fact,  it  is 
shown  in  [ 11 j  that  more  than  one  elastic  and  plastic  region  can  be 
generated  near  the  boundary  even  though  the  stress  state  at  the  boundary 
is  continuously  changing  from  a  lower  yield  surface  to  a  higher  yield 
surface  and  vice  versa.  This  clearly  indicates  the  hidden  difficulties 
in  solving  wave  propagation  of  combined  stresses  by  a  numerical  scheme. 

Ac  know le  dgements 

This  work  was  sponsored  by  the  Ballistic  Research  Laboratories  of 
Aberdeen  Proving  Ground,  Maryland,  through  Contract  DA-04-200-AMC-659(X) 
with  Stanford  University.  Their  support  is  gratefully  acknowledged. 

The  authors  would  also  like  to  thank  Professor  E.H*  Lee  for  valuable 
discussions  during  the  course  of  this  study. 


24 


References 


1.  Bleich,  H.H.  and  Nelson,  I.,  "Plane  Waves  in  an  Elastic -Plastic 
Half -Space  Due  to  Combined  Surface  Pressure  and  Shear,"  J. 

Appl.  Mech.  vol.  33,  1966,  149-158. 

2.  Clifton,  R.J.,  "An  Analysis  of  Combined  Longitudinal  and  Torsional 
Plastic  Waves  in  a  Thin-Walled  Tube,"  Fifth  U.S.  National  Congress 
of  Appl.  Mech.,  June  1966 ,  46>-48o. 

3*  Rakhmatulin,  Kh.  A.,  "On  the  Propagation  of  Elastic -Plastic  Waves 
Owing  to  Combined  Loadings,"  Apnl.  Math.  Mech.  (PMM),  Vol.  22, 
1958,  1079-1088. 

4.  Cristescu,  N-,  "On  the  Propagation  of  Elastic -Plastic  Waves  for 
Combined  Stresses,"  Appl.  Math.  Mech.  (PMM),  Vol.  23,  1959, 

I605-I612. 

5*  Nan,  Ning,  "Elastic -Plastic  Waves  of  Combined  Stresses,"  Tech. 
Rept.  No.  184,  Dept,  of  Applied  Mechanics,  Stanford  University, 
Stanford,  Calif.,  July  1968. 

6.  Ting,  T.C.T.,  "interaction  of  Shock  Waves  Due  to  Combined  Two 
Shear  Loadings,"  Tech.  Rept.  No.  180,  Div.  of  Engineering 
Mechanics,  Stanford  University,  May  1968. 

7.  Kill,  R.,  "The  Mathematical  Theory  of  Plasticity,"  Oxford, 
Clarendon  Press,  1950. 

8.  Courant,  R.  and  Hilbert,  D.,  "Methods  of  Mathematical  Physics," 
Vol.  II,  Interscience  Publishers,  N.J.  1962. 

9*  Jeffrey,  A.  and  Taniuti,  T.,  "Non-Linear  Wave  Propagation  with 
Applications  to  Physics  and  Magneto -hydrodynamic s , "  Academic 
Press,  1964. 

10.  Lipkin,  J.  and  Clifton,  R.J.,  "An  Experimental  Study  of  Combined 
Longitudinal  and  Torsional  Plastic  Waves  in  a  Thin-Walled  Tube," 
Twelth  International  Congress  of  Applied  Mechanics,  held  at 
Stanford  University,  August  1968. 

11.  Ting,  T.C.T.,  "On  the  Initial  Slope  of  Elastic -Plastic  Boundaries 
in  Combined  Longitudinal  and  Torsional  Wave  Propagation,"  Tech. 
Rept.  No.  188,  Dept,  of  Applied  Mechanics,  Stanford 
University,  Sept .  1968 . 


- 


25 


1 


2 


3 


4 


FIGURE  1 


(b) 


.  STRESS  PATHS  FOR  SIMPLE  WAVE  SOLUTIONS 
02-3,  a  -10  (  k/k0  -  1  )Vl,  (3*5 


(°1*°2)/k0 


AH  EXAMPLE  OF  CENTERED  SIMPLE  WAVES  IN 
WHICH  c] ,  c2,  c(,  AND  cs  ALL  APPEAR. 


Security  Classificaiion 


DOCUMENT  CON’ROL  DATA  -  R  8,  D 

(Security  ( latstllc *ti’  n  of  title,  body  ol  hstroc  i  and  indexing  annataiton  must  be  enterec  when^the^overail_re^ort_j_e_<J^'teaUled^ 


\  ORIGINATING  ACTIVITY  (l  >rporate  author) 

Stanford  University 

Stanford,  California  U.S.A.  9^305 

Department  of  Applied  Mechanics 


2*.  REPORT  JECURITY  CUtMf  CATION 


Unclassified 


26.  GROUP 

N.  A. 


3  REPORT  TITLE 


Plane  Waves  Due  to  Combined  Compressive  and  Shear  Stresses  in  a  Half -Space 


4  DESCRIPTIVE  NOTES  (Type  of  report  end  IncUiutvr  date*) 

Interim  Technical 

3-  authorisi  (Firel  name,  middle  initial,  last  name) 


T.C.T.  Ting  and  Ning  Nan 


f  REPORT  DATE 

September,  1968 

7 a.  TOTAL  NO.  OP  PAGES 

25 

7b.  NO.  OP  REPS 

11 

8*.  CONTRACT  OR  GRANT  NO 

|e«.  ORIGINATOR’S  REPORT  NUMBERI3I  S 

DA -04-200 -AMC -65  9 ( X ) 

6.  PROJEC  T  NO 

Technical  Report  No. 

185 

e. 

06.  OTHER  REPORT  NO(3)  (Any  other  numbeta  that  may  be  aaalgned 
thla  report) 

d. 

Contract  Report  No.  12 

10  OIJTRI3UTION  5TATEHEU 


This  docrn*'nt  h'-'”1  -'pnrcvyi  for  public  release 
and  sale;  res  diso^. - jwvica  is  unlimited. 

11  SUPPLEMENTARY  NOTE? 

Hone 

Td  ABSTRACT 

The  plane  wave  propagation  in  a  half -space  due  to  a  uniformly  distributed 
step  load  of  pressure  and  shear  on  the  surface  was  first  studied  by  Bleicl  and 
Nels  n.  The  material  in  the  half-space  was  assumed  to  be  elastic -ideally  plastic. 
In  this  paper,  we  study  the  same  problem  for  a  generil  elastic -plastic  material. 
The  half-space  can  be  initially  prestressed.  The  results  can  be  extended  t.  the 
case  in  which  the  loads  on  the  surface  are  not  necessarily  step  loads,  but  with 
a  restricted  relation  between  the  pressure  and  the  shear  stresses. 


[12.  SPONSORING  MILITARY  ACTIVITY 

tt.  r.  • 

Aberdeen  t---- 
Aberdeen  ivc.c.'  . 


-•1.  Cent: 
.  21, 


DD  ,Fr.,1473 


Sorurilv  ClrtRA-licrflion 


Security  CUsslflceticn 


KCV  WORDS 


Dynamic  Plasticity 

Wave  Propagation  of  Combined  Stresses 
Simple  Wave  Solutions 


